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We consider potential type dynamical systems in finite dimen- 
sions with two meta-stable states. They are subject to two sources of 
perturbation: a slow external periodic perturbation of period T and 
a small Gaussian random perturbation of intensity e, and, therefore, 
are mathematically described as weakly time inhomogeneous diffu- 
sion processes. A system is in stochastic resonance, provided the small 
noisy perturbation is tuned in such a way that its random trajecto- 
ries follow the exterior periodic motion in an optimal fashion, that is, 
for some optimal intensity e(T). The physicists' favorite, measures 
of quality of periodic tuning — and thus stochastic resonance — such 
as spectral power amplification or signal-to-noise ratio, have proven 
to be defective. They are not robust w.r.t. effective model reduc- 
tion, that is, for the passage to a simplified finite state Markov chain 
model reducing the dynamics to a pure jumping between the meta- 
stable states of the original system. An entirely probabilistic notion 
of stochastic resonance based on the transition dynamics between 
the domains of attraction of the meta-stable states — and thus fail- 
ing to suffer from this robustness defect — was proposed before in the 
context of one-dimensional diffusions. It is investigated for higher- 
dimensional systems here, by using extensions and refinements of the 
Freidlin-Wentzell theory of large deviations for time homogeneous 
diffusions. Large deviations principles developed for weakly time in- 
homogeneous diffusions prove to be key tools for a treatment of the 
problem of diffusion exit from a domain and thus for the approach of 
stochastic resonance via transition probabilities between meta-stable 
sets. 
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Introduction. The ubiquitous phenomenon of stochastic resonance has 
been studied by physicists for about 20 years and recently discovered in 
numerous areas of natural sciences. Its investigation took its origin in a 
toy model from climatology, which may serve to explain some of its main 
features. 

To give a qualitative explanation for the almost periodic recurrence of cold 
and warm ages (glacial cycles) in paleoclimatic data, Nicolis [12] and Benzi, 
Sutera and Vulpiani [3] proposed a simple stochastic climate model based on 
an energy balance equation for the averaged global temperature T(t) at time 
t. The balance between averaged absorbed and emitted radiative energies 
leads to a deterministic differential equation for T(t) of the form 



The solar constant Q(t) fluctuates periodically at a very low frequency of 
1CP 5 times per year due to periodic changes of the Earth orbit's eccen- 
tricity (Milankovich cycles), which coincide roughly with the observed fre- 
quency of ice and warm ages. Under reasonable assumptions, for frozen q, 
the nonlinear function b(q,T) describes the force associated with a double 
well potential possessing two stable temperature states which represent cold 
and warm ages. As Q varies periodically, these states become meta-stable 
and are moved periodically by Q(t). Most importantly, transitions between 
these states are impossible. Only the addition of a stochastic forcing allows 
for spontaneous transitions between the meta-stable climate states, thus ex- 
plaining roughly transition mechanisms leading to glacial cycles. 

In general, trajectories of the solutions of differential equations of this 
type, subject to two independent sources of perturbation, an exterior pe- 
riodic one of period T, and a random one of intensity e, say, will exhibit 
some kind of randomly periodic behavior, reacting to the periodic input 
forcing and eventually amplifying it. The problem of optimal tuning at large 
periods T consists in finding a noise amplitude e(T) (the resonance point) 
which supports this amplification effect in a best possible way. During the 
last 20 years, various concepts of measuring the quality of periodic tuning 
to provide a criterion for optimality have been discussed and proposed in 
many applications from a variety of branches of natural sciences (see [7] for 
an overview). Its mathematical treatment started only very recently, and 
criteria for finding an optimal tuning are still under discussion. 

The first approach toward a mathematically precise understanding of 
stochastic resonance was done by Preidlin [5] . Using large deviations theory, 
he explains basic periodicity properties of the trajectories in the large period 
(small noise) limit by the effect of deterministic quasi-periodic motion, but 
fails to account for optimal tuning. The most prominent quality measures 
for periodic tuning from the physics literature, the signal-to-noise ratio and 




T(t) = b(Q(t),T(t)). 
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the spectral power amplification coefficient (SPA), were investigated in a 
mathematically precise way in Pavlyukevich's thesis [13], and seen to have a 
serious drawback. Due to the high complexity of original systems, when cal- 
culating the optimal noise intensity, physicists usually pass to the effective 
dynamics of some kind of simple caricature of the system reducing the dif- 
fusion dynamics to the pure inter well motion (see, e.g., [11]). The reduced 
dynamics are represented by a continuous time two state Markov chain. 
Surprisingly, due to the importance of small intra well fluctuations, the tun- 
ing and resonance pattern of the Markov chain model may differ essentially 
from the resonance picture of the diffusion. It was this lack of robustness 
against model reduction which motivated Herrmann and Imkeller [9] to look 
for different measures of quality of periodic tuning for diffusion trajectories, 
retaining only the rough inter well motion of the diffusion. The measure 
they treat in the setting of one-dimensional diffusion processes subject to 
periodic forcing of small frequency is related to the transition probability 
during a fixed time window of exponential length, the position of which is 
tracked by a parameter of period length in which maximization is performed 
to account for optimal tuning. 

The subject of the present paper is to continue our previous work in 
the general setting of finite- dimensional diffusion processes. Our approach 
of stochastic resonance thereby is based on the same robust probabilistic 
notion of periodic tuning. This extension is by no means obvious, since 
the multidimensional problem requires entirely new methods. We recall at 
this point that in [9] methods of investigation of stochastic tuning were 
heavily based on comparison arguments which are not an appropriate tool 
from dimension 2 on. Time inhomogeneous diffusion processes, such as the 
ones under consideration, were compared to piecewise homogeneous diffusion 
processes by freezing the potential's time dependence on small intervals. We 
study a dynamical system in c£-dimensional Euclidean space perturbed by a 
c?-dimensional Brownian motion W, that is, we consider the solution of the 
stochastic differential equation 



One of the system's important features is that its inhomogeneity is weak 
in the sense that the drift depends on time only through a re-scaling by 
the time parameter T = T e which will be assumed to be exponentially large 
in s. This corresponds to the situation in [9] and is motivated by the well- 
known Kramers-Eyring law which was mathematically underpinned by the 
Freidlin-Wentzell theory of large deviations [6] . The law roughly states that 
the expected time it takes for a homogeneous diffusion to leave a local attrac- 
tor, for example, across a potential wall of height |, is given to exponential 
order by T e = exp(|). Hence, only in exponentially large scales of the form 



(0.2) 




t > 0. 
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T £ = exp(^) we can expect to see effects of transitions between different 
domains of attraction, b is assumed to be one-periodic w.r.t. time. The de- 
terministic system £ t = b(s,£t) with frozen time parameter s is supposed to 
have two domains of attraction that do not depend on s > 0. In the "classi- 
cal" case of a drift derived from a potential, b(t,x) = —V x U(t,x) for some 
potential function U, equation (0.2) describes the motion of a Brownian 
particle in a <i-dimensional time inhomogeneous double-well potential. 

Since our stochastic resonance criterion is based on transition times be- 
tween the two meta-stable sets of the system, our analysis relies on a suitable 
notion of transition or exit time. The Kramers-Eyring formula suggests to 
consider the parameter \x from T £ = exp(^) as a natural measure of scale. 
Therefore, if at time s the system needs energy e(s) to leave some meta- 
stable set, an exit from that set should occur at time 

a M = inf{t > : e(t) < fx} 

in the diffusion's natural time scale. If a], are the transition times for the 
two domains of attraction numbered i = ±1, we look at the probabilities of 
transitions between them within a time window [(ct* — h)T £ , (a\ + h)T £ ] for 
small h > 0. Assume for this purpose that the two corresponding meta-stable 
points are given by Xi,i = ±1, and denote by t~ % the random time at which 
the diffusion reaches the ^-neighborhood B g {x^i) of x_j. Then we use the 
following quantity to measure the quality of periodic tuning: 

M(e^) = mm sup F x (t~ 1 G [(a* - h)T, (a* + h)T}), 
t=±1 xeB e {xi) 

the minimum being taken in order to account for transitions back and forth. 
In order to exclude trivial or chaotic transition behavior, the scale parameter 
/i has to be restricted to an interval Ir of reasonable values which we call 
resonance interval. With this measure of quality, the stochastic resonance 
point may be determined as follows. We fix the window width parameter 
h > and maximize M(e,n) in \x asymptotically as e — > 0, or, equivalently, 
as T £ = exp j — ► oo. If the asymptotic maximizer exists and is reached for 
the time scale /J-o(h), we call the eventually existing limit /j,q = lim^o [J>o(h) 
resonance point. This is the optimal tuning w.r.t. our quality measure, that 
is, the best asymptotic relation between the noise amplitude e and the period 
T £ . 

To calculate Ho(h) for fixed positive h, we use large deviations techniques. 
In fact, our main result contains a formula which states that 

lim elog{l — -M(e, fj,)} = max{/i — ei(a l — h)}. 

e^O i=±l ^ 

We show that this asymptotic relation holds uniformly w.r.t. (J, on compact 
subsets of Ir, a fact which enables us to perform a maximization and find 
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fj>o(h). The techniques needed to prove our main result feature extensions 
and refinements of the fundamental large deviations theory for time ho- 
mogeneous diffusions by Freidlin-Wentzell [6]. We prove a large deviations 
principle for the inhomogeneous diffusion (0.2) and strengthen this result 
to get uniformity in system parameters. Similarly to the time homogeneous 
case, where large deviations theory is applied to the problem of diffusion 
exit culminating in a mathematically rigorous proof of the Kramers-Eyring 
law, we study the problem of diffusion exit from a domain which is carefully 
chosen in order to allow for a detailed analysis of transition times. The main 
idea behind our analysis is that the natural time scale is so large that re- 
scaling in these units essentially leads to an asymptotic freezing of the time 
inhomogeneity, which has to be carefully studied, to hook up to the theory 
of large deviations of time homogeneous diffusions. 

The material in the paper is organized as follows. Section 1 is devoted 
to the careful extension of large deviations theory to diffusions with slow 
time inhomogeneity. The most useful result for the subsequent analysis of 
exit times is Proposition 1.8, with a large deviations principle for slowly 
time dependent diffusions, uniform with respect to a system parameter. In 
Section 2 upper and lower bounds for the asymptotic exponential exit rate 
from domains of attraction for slowly time dependent diffusions are derived. 
The main result Theorem 2.3 combines them. Section 3 is concerned with 
developing the resonance criterion and computing the resonance point from 
the results of the preceding section. 

1. Large deviations for diffusion processes. Let us now consider dynam- 
ical systems driven by slowly time dependent vector fields, perturbed by 
Gaussian noise of small intensity. We shall be interested in their large devi- 
ation behavior. Due to the slow time inhomogeneity, the task we face is not 
covered by the classical theory presented in [4] and [6]. For this reason, we 
shall have to extend the theory of large deviations for randomly perturbed 
dynamical systems developed by Freidlin and Wentzell [6] to drift terms de- 
pending in a weak form to be made precise below on the time parameter. 
Before doing so in the second subsection, we shall recall the classical results 
on time homogeneous diffusions in the following brief overview. 

1.1. The time homogeneous case: classical results. For a more detailed 
account of the following well-known theory, see [4] or [6]. 

We consider the family of Revalued processes X s , e > 0, defined by 

(1.1) dXf = b(Xf)dt + ^dW t , X^ = x G R d , 

on a fixed time interval [0, T], where b is Lipschitz continuous and W is a 
d-dimensional Brownian motion. This family of diffusion processes satisfies 
in the small noise limit, that is, as e — > 0, a large deviations principle (LDP) 
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in the space Cqt '■= C([0,T],M. d ) equipped with the topology of uniform 
convergence induced by the metric poT(p,ip) ■= sup 0<t<T \\tpt — ipt\\, ip, ip £ 
Cqt- The rate function or action functional is given by Iq£ : Cqt [0, +00], 



\ [ 7 \\<pt - b((f t )f dt, 

(1.2) I&(<p) = { Jo - f u 1,1 y j 

UJ 11 ip is absolutely continuous and ip§ = xq , 

+00, otherwise. 

Moreover, I^j, is a good rate function, that is, it has compact level sets. 
The LDP for this family of processes is mainly obtained as an application 
of the contraction principle to the LDP for the processes y/sW, e > 0. More 
precisely, in the language of Freidlin and Wentzell, the functional Jqj, is the 
normalized action functional corresponding to the normalizing coefficient p 
In the sequel we will not consider scalings other than this one. We have 
Iqt^) < 00 if an d only if ip belongs to the Cameron-Martin space of ab- 
solutely continuous functions with square integrable derivatives starting at 
xo, that is, 

<peH$:=[f:[0,T}^m d \f(t)=x + J\(s)ds for some g E L 2 ([0,T])j. 

We omit the superscript xq whenever there is no confusion about the initial 
condition we are referring to. 

Observe that Iot(v) = means that (p (up to time T) is a solution of the 
deterministic equation 

(1.3) i=Ko, 

so Iot(<p) is essentially the L 2 -deviation of ip from the deterministic solution 
£. The cost function V of X s , defined by 

V{x, y, t) = mf{I ot (ip) :peC ot ,po = x,p t = y}, 

takes into account all continuous paths connecting x, y S W 1 in a fixed time 
interval of length t, and the quasi-potential 

V(x,y) = inf V(x,y,t) 

describes the cost of X £ going from x to y eventually. In the potential case, 
V agrees up to a constant with the potential energy to spend in order to 
pass from x to y in the potential landscape, hence, the term quasi-potential. 

The classical LDP due to Freidlin and Wentzell requires the usual global 
Lipschitz and linear growth conditions from the standard existence and 
uniqueness results for SDE. In our setting the coefficients will (in general) 
not be globally Lipschitz since the drift is given by a potential gradient. An 
extension to locally Lipschitz and e-dependent drift terms was provided by 
Azencott [1]. The following proposition is a special case of [1], Chapter III, 
Theorem 2.13. See also [2], Theorem 2.1. 
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Proposition 1.1. Assume that the equation (1.1) has a unique strong 
solution that never explodes and that the drift is locally Lipschitz. Then X s 
satisfies on any time interval [0, T] a large deviations principle with good 
rate function Iqt- Furthermore, the LDP for X s holds uniformly w.r.t. the 
initial condition of the diffusion. More precisely, if F y (X £ G •) denotes the 
law of the diffusion X s starting in y G M rf and K C M. d is compact, we have, 
for any closed F C Cot, 

(1.4) limsupelogsupP^pT G F) < - inf inf lZr{<P) 

e^o y eK yeKipeF 

and for any open G C Cot, 

(1.5) liminf elog inf F y (X £ G G) > - sup inf iKrfa). 

Remark 1.2. (i) A sufficient condition for the existence of a nonex- 
ploding and unique strong solution is a locally Lipschitz drift term b which 
satisfies 

(1.6) (x,6(x))<7(l + ||x|| 2 ) forallxGK^ 

for some constant 7 > (see [17], Theorem 10.2.2). This still rather weak 
condition is obviously satisfied if (x,b(x)) <0 for large enough x, which 
means that b contains a component that pulls X back to the origin. 

(ii) A strengthening of condition (1.6) ensuring superlinear growth will be 
used in subsequent sections. In that case, the laws of {X £ ) are exponentially 
tight, and Iqt is a good rate function. Recall that the laws of {X s ) are 
exponentially tight if there exist some Rq > and a positive function ip 
satisfying lim^^oo (p(x) = +00 such that 

(1.7) limsupelogP(crfj < T) < -<p(R) for all R > R . 

Here a e R denotes the first time that X £ exits from Br(0). 

1.2. General results on weakly time inhomogeneous diffusions. Let us 
now come to inhomogeneous diffusions with slowly time dependent drift co- 
efficients. For our understanding of stochastic resonance effects of dynamical 
systems with slow time dependence, we have to adopt the large deviations 
results of the previous subsection to diffusions moving in potential land- 
scapes with different valleys slowly and periodically changing their depths 
and positions. In this subsection we shall extend the large deviations results 
of Freidlin and Wentzell to time inhomogeneous diffusions which are almost 
homogeneous in the small noise limit, so that, in fact, we are able to com- 
pare to the large deviation principle for time homogeneous diffusions. The 
result we present in this subsection is not strong enough for the treatment 
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of stochastic resonance (one needs uniformity in some of the system param- 
eters), but it most clearly exhibits the idea of the approach, which is why 
we state it here. Consider the family X £ , e > 0, of solutions of the SDE 

(1.8) dX\ = b £ (t, XI) dt + V^dWt, t>0, X £ = x € M. d . 

We assume that (1.8) has a global strong solution for all e > 0. Our main 
large deviations result for diffusions for which time inhomogeneity fades out 
in the small noise limit is summarized in the following proposition. The e- 
dependence of the drift term was assumed in the same way in [1], Chapter 
III, Theorem 2.13 and [2], Theorem 2.1. See also [14]. 

Proposition 1.3 (Large deviations principle). Assume that the drift of 
the SDE (1.8) satisfies 

(1.9) lim b £ {t,x) = b(x) 

for all t>0, uniformly w.r.t. x on compact subsets ofM. d , for some locally 
Lipschitz function b:M. d —>U. d . Assume that the time homogeneous diffusion 
Y £ associated to the limiting drift b [i.e., the solution of (1.1) with the same 
initial condition xo] does not explode. 

Then (X £ ) satisfies a large deviations principle on any finite time interval 
[0,T] with good rate function Iqt given by (1.2). 

Proof. For notational convenience, we drop the e-dependence of X and 
Y . We shall prove that X and Y are exponentially equivalent, that is, for 
any 5 > 0, we have 

(1.10) limsupelogF( PoT (X,Y) > 5) = -co. 

e^O 

In order to verify this, fix some 5 > 0, and observe that 

||^-^||< f t \\b £ (u,X u )-b(X u )\\du+ f \\b{X u )-b(Y u )\\du. 
Jo Jo 

For R > 0, let tr := inf{t > 0:X t ^ Br(xo)}, let ?r be defined similarly 
with X replaced by Y, and aR := tr Atr. The local Lipschitz continuity of 
b implies the existence of some constant Kr(xq) such that \\b(x) — b(y)\\ < 
Kr(xo)\\x — y\\ for x, y G Br(xq). An application of Gronwall's lemma yields 

P0T{X,Y)<e KR ^ T [ T \\b £ (u,X u )-b(X u )\\du on{a R >T}. 
Jo 

Due to uniform convergence, for any rj > 0, we can find some Eq > s.t. 
sup \\b £ (t,x) -b(x)\\ <7] for te [0,T],e<£ - 

x£B R (x ) 
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This implies 

(1.11) p T(.X,Y)<r]Te KR(xo)T for e < e on {a R > T}. 

By choosing rj small enough s.t. pot(X,Y) < 5/2 on {a R > T} [i.e., X and 
Y are very close before they exit from B R (xo)], we see that, for e < eq, 

F(p 0T (X, Y)>6)< F(t r < T) + F(f R < T). 

Since X and Y are close within the ball B r (xq), we deduce that if X escapes 
from B r (xq) before time T, then Y must at least escape from B r /2{xq) 
before time T (if R > 5). So we have ¥(p 0T (X,Y) > 5) < P(f fl/2 < T) for 
e < £q. Hence, the LDP for Y gives 

limsupelogP(o T(^,^) >^) < - inf V(x ,y,t). 

e^O 0<t<T\\y-x \\>R/2 

Sending R — > oo yields the desired result (see Theorem 4.2.13 in [4]). □ 

It is easy to see that the uniformity w.r.t. the diffusion's initial condition 
also holds for the weakly inhomogeneous process X s of this proposition. One 
only has to carry over Proposition 5.6.14 in [4], which is easily done using 
some Gronwall argument. Then the proof of the uniformity is the same as 
in the homogeneous case (see [4], Corollary 5.6.15). We omit the details. 

1.3. Weak inhomogeneity through slow periodic variation. In this sub- 
section we shall deal with some particular diffusions for which the drift 
term is subject to a very slow periodic time inhomogeneity. More precisely, 
we shall be concerned with solutions of the following stochastic differential 
equation taking their values in <i-dimensional Euclidean space, driven by a 
(i-dimensional Brownian motion W of intensity e: 

(1.12) dX £ t =b^,X^dt + ^dW t , t>0,X = x eM. d . 

Here T £ is a time scale parameter which tends to infinity as e — > 0. In the 
subsequent sections, we shall assume that T 6 is exponentially large, in fact, 

(1.13) T £ = exp^ with / u>0. 

The drift b(t,x) of (1.12) is a time-periodic function of period one. Concern- 
ing its regularity properties, we suppose it to be locally Lipschitz in both 
variables, that is, for R > 0, i£ ]R d , there are constants K R (x) and k r (x) 
such that 

(1-14) \\b(t,y 1 )-b(t,y 2 )\\<K R (x)\\y 1 -y 2 \\, 
(1.15) \\b(t,y)-b(s,y)\\<K R (x)\t-s\ 
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for all y, 2/1, 2/2 G Br(x) and s,t>0. Furthermore, we shall assume that the 
drift term forces the diffusion to stay in compact sets for long times in order 
to get sufficiently "small" level sets. We suppose that there are constants r], 
Ro > such that 

(1.16) (x,b(t,x)) < —7]\\x\\ 

for t > and ||x|| > Ro. This condition is stronger than (1.6), so the existence 
of a unique strong and nonexploding solution is again guaranteed. Moreover, 
this growth condition implies the exponential tightness of the diffusion (see 
Proposition 1.4 for the precise asymptotics) . 

1.3.1. Boundedness of the diffusion. The aim of this subsection is to 
exploit the consequences of the growth condition (1.16). In fact, it implies 
that the diffusion (1.12) cannot leave compact sets in the small noise limit. 
For positive e, it stays for a long time in bounded domains. In the following 
proposition we shall make precise how the law of the exit time from bounded 
domains depends on e. The arguments are borrowed from the framework of 
self-attracting diffusions; see [15] or [10]. 

For R > 0, let a £ R := inf{i > : \\Xf || > R} denote the first exit time from 
the ball B R (0). 

Proposition 1.4. Let 5 > 0, and let r:(0,5) — > (0,oo) be a function 
satisfying lim £ ^o ffey = 0. There exist R\,£\ > and C > such that, for 

R > R\ , £ < £\, 

(1.17) F x (a £ R < r(e)) < CV lM e -^A f or || x || < 

Proof. For convenience of notation, we suppress the superscript e in 
X £ , a £ R and so on. Choose a C 2 -function h : M. d — > R s.t. h(x) = \\x\\ for ||x|| > 
Rq and h(x) < Rq for ||x|| < Rq, where Rq is the constant given in the growth 
condition (1.16). By Ito's formula, we have 

h(X t ) = h(x) + f Vh(X s ) dW s 
Jo 

+ J*(vh,b(j^, •) \(X S ) ds + e -J^ Ah(X s ) ds. 

Let '■= Jo ||V/i(X s )|| 2 ds, that is, £t is the quadratic variation of the con- 
tinuous local martingale M t := J Q Vh(X s ) dW s ,t > 0. Since Vh(x) = for 
\\x\\ > Rq, we have d£t = dt on > Ro}. Now we introduce an auxiliary 

process Z which shall serve to control ||X||. 
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Let < fj < rj. According to Skorokhod's lemma (see [16]), there is a 
unique pair of continuous adapted processes (Z, L) such that L is an in- 
creasing process (of finite variation) which increases only at times t for which 
Z% = Rq , and Z > Rq, which satisfies the equation 

Z:=R V \\x\\ + JeM -fj£ + L. 

We will prove that 

(1.18) \\X t \\<Z t a.s. for all t > 0. 
For that purpose, choose / € C 2 (K) such that 

f(x) > and f'(x) > for all x > 0, 
f(x) = for all x < 0. 

According to Ito's formula, for t > 0, 

f(h(X t ) - Z t ) = f(h(x) - \\x\\ V Rq) + / f'(h(X s ) - Z s ) d(h(X) - Z) s 

Jo 

+ \ f f"(h(X s ) - Z s ) d(h(X) - Z) s . 
Jo 

By definition of h and Z, we have h(Xt) < Zt on < Ro}, so {h(Xt) > 

Zt} = {\\Xt\\ > Zt}. Moreover, by definition, h(X) — Z is a finite variation 
process. Hence, the expression 

- f f{\\X s \\-Z s )dL s 
Jo 

is an upper bound of f(h(X t ) — Z t ). Furthermore, Ah(x) = for ||x|| > Rq, 
which by (1.16) implies 

j±j ^X s ,b(^,X^y £ -Ah(X s )+f}< £ -^^+fj-r, on{||X s ||>Z,}. 

The latter expression is negative if e is small enough, so we can find some 
£o > such that — Z t ) < for e < £q. This implies \\X t \\ < Z t a.s. by 

the definition of /, and (1.18) is established. 

We therefore can bound the exit probability of X by that of Z. If Q 
denotes the law of the process Z, we see that, for any a > 0, 

(1.19) F x (a R < r(e)) < Q(a R < r(e)) < e^^Egfe"^]. 

In order to find a bound on the right-hand side of (1.19), let us define 
K := sup, W j <R ||V/i(x)|| 2 . Then we have & < Kt for all t > 0. Note that 
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w.l.o.g. h can be chosen so that K < 2Rq. Now observe that, by Ito's formula, 
for any ip£C 2 (R), 

d(^(Z t )e-^ K ^) = y/Ztf/iZje-WW* d M t + p'(Z t ) e -^ K ^ dL t 

+ e-(«/^{| v "(Z t )-^(Z t ) - %<p(Z t )}d£ t . 
Now let R> Rq. If we choose <p such that 

E O. 

-tp (y) - fj(p(y) - —ip(y) = for y £ [R ,R], 

<p'(R ) = 0, <p(R) = 1, 

then (p(Z t )e~( a / K ^ t is a local martingale which is bounded up to time a R . 
Hence, we are allowed to apply the stopping theorem to obtain 

(1.20) <p(R V \\x\\) = E Q [<p(Z aR )e-W K ^n] = E Q [e^ a / K ^R]. 

Hence, & R < Ka R , which implies E Q [e- (a/i ^ CT fl] > E Q [e- aaR ], and we de- 
duce from (1.19) that 

(1.21) F x (a R < r(e)) < e ar ^E Q [e^ K ^R] < e ar ^p(R V ||z||). 

Solving the differential equation for tp yields 

_\~ e \ + {x-R ) _|_ \+ e \~{x-Ro) 

with A ^ j^» , Hence, 

, , (A+-A-)e A+ ( a; - Ko ) 

^ X) - {-x-^nR-Ro) ■ 

Taking a = r(e)~ 1 in (1.21), we obtain 

W x (a R < r(e)) < exp(l)^( J R V ||x||) < A+ ~ j exp{l + A+(i? V - #)}. 

— A 

It is obvious that exp{A+(i? V ||x|| - J?)} < exp{-^} for R > 2(\\x\\ V Rq), 
so it remains to comment on the prefactor. We have 

A+-A- 2y/fj* + 2(a/K)e 4(f? 2 + 2e/(Kr{e))) 
-A" ~7^ 2 + 2(a/^)e^- 2e/(iO( e )) 

Since lim E _»o = 0, the latter expression behaves like 2fj 2 K^-L as e — > 0. 
Putting these estimates together yields the claimed asymptotic bound with 
fj instead of rj. Letting fj — > rj completes the proof. □ 
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Remark 1.5. The proof of Proposition 1.4 shows a lot more. The cru- 
cial inequality (1.21) contains a bound which is independent of X £ , since ip 
is defined by means of h, e, fj and Rq only. Thus, we have shown that the 
bound (1.17) holds for all diffusions satisfying the growth condition (1.16), 
that is, e\ and R\ are independent of X s . In particular, (1.17) holds uni- 
formly w.r.t. fx. 

1.3.2. Properties of the quasi-potential. Taking large period limits in the 
subsequently derived large deviations results for our diffusions with slow 
periodic variation will require to freeze the time parameter in the drift term. 
The corresponding rate functions are given a separate treatment in this 
subsection. We shall briefly discuss their regularity properties. This will 
be of central importance for the estimation of exit rates in Section 2. For 
s > 0, T > 0, we consider 



As in the first section, we need associated cost functions. For s > 0, x, y G M. d , 
they are given by 



to be at the point y at time t when starting at x. The corresponding quasi- 
potential 



describes the cost for the frozen system to go from x to y eventually. Let 
us note that since the drift b is locally Lipschitz in the time variable, the 
family of action functionals Iq T is continuous w.r.t. the parameter s, and the 
corresponding cost functions and pseudo-potentials inherit this continuity 
property. Let us recall some further useful properties of the quasi-potentials 
and their underlying cost and rate functions. The following properties are 
immediate. 

Lemma 1.6. For any x,y,z G IR d and s,t,u> 0, we have the following: 



(1.22) 




if if is absolutely continuous 



(1.23) V s (x,y, t) = inf{/ot(y) : <p G C ot ,<Po = x,tp t = y}. 

V s (x,y,t) is the cost of forcing the frozen system 

dY t £ = b(s,Y t £ )dt + ^dW t , t>0, 



(1.24) 



V s (x,y)=mfV s (x,y,t) 



(a) V s (x,y,t + u) <V s (x,z,t) + V s (z,y,u), 

(b) (s,y) i— > V s (x,y,t) is continuous on x IR'', 

(c) inf|i„|i>^ V s (x, y, t) — ► oo uniformly w.r.t. s > 0. 
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The following lemma establishes the local Lipschitz continuity of the 
quasi-potential w.r.t. the state variables, uniformly w.r.t. the parameter s. 

Lemma 1.7. For any compact subset K ofM. d , there exists Tk > such 
that 

supV s (x,y) < r^dist(a;, y) 

s>0 

for all x,y € K . 

Proof. Let x and y belong to K. There exists some radius R > such 
that K C B R (0). Set T = dist(x,y). We construct a path ip £ Cot by setting 
(ft = x + di ^xy) t ^ or * e [0>^T Since 6(s, •) is locally Lipschitz uniformly 
w.r.t. s > 0, we obtain an upper bound for the energy of 93: 

1 /" T 

Io T (p) < -sup \\(p t - b(u,ip t )\\ 2 dt 

* u>0 JO 

2 Jo \dist(x,y) <«<i / 
<\ / T (l + Kfl (0) + ||6(0,^)||) 2 ^ 

< \ f T (l + K R (0) + K R (0)y t \\ + ||6(0,0)||) 2 (it 
2 Jo 

< |(1 + Kfl(0) + RK R (0) + ||6(0,0)||) 2 . 

For T K := ±(1 + k r (0) + RK R (0) + ||6(0,0)||) 2 and by the definition of T, 
we obtain 

supV s (x,y) < sup Iq T (ip) < r K dist(x,y). 

s>0 s>0 u 

1.3.3. Large deviations. We shall now specialize the general large devia- 
tions results of the previous subsection to the family X s , e > 0, of solutions 
of (1.12). At the same time they will be strengthened, to obtain uniformity 
w.r.t. some of the system's parameters: the scale parameter //, the starting 
time and the initial condition. 

It is an immediate consequence of Proposition 1.3 that the solution of (1.12) 
satisfies a large deviations principle with rate function I® T , that is, the 
rate function is the same as that of a homogeneous diffusion governed by 
the frozen drift 6(0,-). In order to see this, one only has to mention that 
lim £ _ > o&(^->^) =6(0, sc) locally uniformly w.r.t. x due to the Lipschitz as- 
sumptions on b. 
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But this result is not strong enough. We also need some uniformity w.r.t. 
the starting times of the diffusions we consider. Our large deviations state- 
ments derived so far rely on comparison arguments which yield exponen- 
tial equivalence with time homogeneous diffusions for which an LDP is well 
known from the classical theory of Freidlin and Wentzell. In order to achieve 
uniform large deviations estimates, we have to refine this technique, to de- 
rive a large deviations principle for our family of diffusions (1.12), which 
is uniform with respect to both the starting time and the scale parameter. 
This will be our main tool for estimating the asymptotics of exit time laws 
in the subsequent section. 

The diffusion (1.12) is a time inhomogeneous Markov process. The solu- 
tion starting at time r > with initial condition x £ M. d has the same law as 
the solution X r ' x of the SDE 

(1.25) dXl' x = b{^,X r t ^ dt + y/idWt, t > 0, X r ' x = xeR d . 

We denote its law by IP^rO), assume from now on that T e = exp j for some 
fi > 0, and fix T > 0. 

Proposition 1.8. Let K cR d be a compact set and V C (0,oo). For 
/* £ V, r £ [0, 1] and /3>0, let S r, P(e,[i) be a neighborhood of rT e such that 

diam(S^( £ ,^)) ^ 
hmsup sup — < p. 

e^O ^ev,re[o,i] T £ 

Then for any closed F C Cqt, there exists 5 = 5(F) such that 

limsupelog sup F yu (X e eF) < — inf inf Iq T (<^), 

where 70 = Jo(F) = (35(F) and F 10 is the closed ^-neighborhood of F . For 
any open G C Cqt, there exists 5 = 5(G) and (3q = (3q(G) such that, if [3 < po, 

liminfelog inf F vu (X e £G) > — sup inf In T (tp), 

where 70 = 70 (GO = (35(G) and G 10 is the complement of (G c ) 70 . 
These bounds hold uniformly w.r.t. r. 

Remark 1.9. The upper bound means that, for any t? > 0, we can find 
£0 > s.t. for e < £0, we have 

elog sup P y , M (X £ GF)<-inf inf VM+i 

y eK, f ,eV,ueS r -i 3 (e )A ») y &K <P&F~m, V0 =y 

The uniformity in the statement means that £0 can be chosen independently 
of r. A similar statement holds for the lower bound. 
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Observe that the expression for the blowup-factor 70 (F) depends on the 
set F only through 5(F) which is independent of (3, and that 70 (F) — > as 
(3 — > for all F. In particular, if (3 is equal to zero, we recover the classical 
bound of the uniform LDP. 

Proof of Proposition 1.8. For y G R d and r G [0, 1], let Y r ' y be the 
solution of the homogeneous SDE 

dY t r ' y = b(r, Yp y ) dt + y/ZdW u t > 0, Y^ y = y. 

Let W C [0, 1] and r G W. For r G W, u G 5 r '^(e, fJ,),fMG V and i? > 0, let 
:= inf{t > : £ S fl (0)}, and let f£?' y be defined similarly with X u ^ 
replaced by Y r ^ , and <r]f ' r ° := rjf A f^. 

As a consequence of Gronwall's lemma, we see just as in the proof of 
Proposition 1.3 that, for r, tq G [0, l],u G S r,l3 (e,fi), 



p 0T (X u > y ,Y ro > y )<e KRi0)T [ 

Jo 

on {(Jr V ' T0 > T}. This implies 



dt 



P ot(X^,Y^) < . fl (0)re^(^( diam % ( ^ ))+T + |r - r 

on {<7^' y ' r ° > T}. Due to our assumption, the last expression is bounded by 

(1.26) 0i = /3i(W)=p o (W)M(R) ase^O, 
where 

(1.27) Pq(W) :=(3+ sup |r - r | and M(i?) := TK R (0)e KR{0)T . 

r-GW 

Upper bound. Fix some closed set F C Cqt- For all 7 > 0, we have 
F(X u ' y G F) < P(y r °' J/ g F"<) + F( PoT (X u > y ,Y r °> y ) > 7). 
This yields 

limsupelog sup Fy )U (X £ G F) 

yGiC,/iGV,rGW,MeS r >' 3 (e,M) 

(1.28) < lim sup e log max ( sup F(Y r ° ' y G F 7 ), 

sup P(poT(^,r ro ^)>7)). 

1/ei ; !',AtGV,reW,uGS'-.' 3 (£,^) J 
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Now we wish to find some 7 such that the asymptotics of the maximum is 
determined by the left term sup yGK ¥(Y r °' v G F 1 ). In that case the uniform 
LDP for Y will give us the bound 

limsupelog sup P |)i „(I E e F) 

e^O yeK,n€V,reW,ueS r >P(e tt i) 

(1.29) 

< - inf inf E&Up). 

Unfortunately, such a 7 will depend on F. In order to see that it exists and 
can be chosen as claimed in the statement, we define 

G(R,e):= sup F(t% v <T)+ sup P(t% v <T). 

re[0,l], y eK,ij.eV,ueS r ^(s >t j,) re[0,l],yeK 

By Proposition 1.4 and Remark 1.5, we have limsup^Q elog 0(i?, e) < — r/R 
for all R>R\. Hence, we may find R> R\ such that 

lim sup e log (ii, e) < — sup inf inf Iq T (ip). 

e^o re[o,i] y^ K <peF,tpo=y 

Fix this R, let 5(F) = M(R), and note that 5(F) is independent of (3 and W. 
By (1.26), for any 7 > /?i(W) = /3 (W)5(F), we can find e > such that, 
for e < So , 

(1.30) sup W(p QT (X^,Y r ™) > 7) < @(R,s). 

Now we exploit the definition of 5(F), which ensures that the maximum 
in (1.28) is given by the left term. Indeed, for 7 > /?i(W), we have 

limsupelog sup F(p 0T (X u ' v , Y r °' y ) > 7) 

< limsupelogG(i?, e) < — sup inf inf ^ot(v) 

e^O re[0,l] W&K <f>£F,<p =y 

< — inf inf I^(ip), 

yeK V >eFt m =y UJ 

which implies (1.29). The particular choice W = {ro} yields this bound for 
all 7 > 70(F) = (35(F) given in the statement and proves the claimed bound. 
By taking the limit 7 — > 70(F), we get the asserted upper bound since / is 
a good rate function (see [4], Lemma 4.1.6). 

It remains to prove the uniformity w.r.t. r. For that purpose, fix ■& > 0, and 
for ro G [0, 1], let W ro be a neighborhood of ro. Choose (ft* G F 70 [here 70 = 
70(F) = (35(F)] starting at some point yo G K with ini y& K ^ipeF^o,ip = y Iqt^) — 
Iqj,((()*) — t3/8. By Lemma 4.1.6 in [4] and the continuity of Iq T (<P*), 
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we can assume W ro to be small enough such that, for r G W ro , denoting 
7* = /3o(W ro )<S(n 

inf inf I r °(tp)> ml inf I r ° Up) - 0/8 
>I r 0T (<f>*)-#/2 

Due to compactness, we can choose finitely many points n, --^rjv such that 
their corresponding neighborhoods cover [0,1]. Denote 7* := /3o(W rjl )<5(F). 
For each 1 < n < N, there exists some e n > such that, for e < e n and 
rG VW r „, 

elog sup W y>u (X E G F) < - inf inf i^^) + 1 

y&K,iieV,u&S r (e,/i) yeK ip&F^-n, ,tp =y * 

<-inf inf J£ T (<p)+i?. 

yeK <peFlo ,tp =y 

Hence, for e < mini< n <Are n , the preceding inequality holds for all r £ [0, 1]. 

Lower bound. Let G C Cqt be an open set. Consider the increasing func- 
tion 

/(Z):=-sup inf / o r f(0), 
V y£K<l>eG l :<l> =y 

let Iq = inf{/ > : f(l) = +00}, and recall that r\ is the constant introduced 
in the growth condition (1.16) for the drift. Assume first that Iq < 00 (this 
is guaranteed if G is bounded), and set 

H:=/((io-A)(W))v|)vl2i and 7 := /3 (W)M(R), 

where R± is given by Proposition 1.4. Then 

p(yro,j/ G £ 7 ) < P (X"'f G G) +P(poT(^ r0 ' y ,^ U,S/ ) > 7). 

By the uniform LDP for Y ro,y and (1.30), we conclude that 

-17/(7) = -sup mf /^(^ <liminfelog inf P(r° l!, eG T ) 
y tzK<p&3"fm=y £ ~*° yt K 

< maxlliminfelog inf F(X u,y G G), 

{ e^O r€W,y£K,ij,£V,u€S r (E,n) 

limsupelog sup ¥(p 0T (Y r °' y ,X u ' y ) > 7)) 

e^O r£W,y£K,ij,€V,u£S r (E,Li) > 
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< maxjliminfelog inf F(X u,y G G?), 

t £^0 reW,yGK,iJ,GV,u€S r (e,fi) 

lim sup £ log Q(R, e) 

£-►0 

Since / is increasing and R> R±, we obtain by Proposition 1.4 
-r//(7 + /3o(W)) 

< -vfd) 

< max (lim inf e log inf ¥(X u ' y G G), -vr\- 

{ r£W,y£K,/j,£V,u£S r (e,n) J 

Now we have to compare f("f) and R in order to see when the maximum is 
given by the left term. 

If /(7) > R, then 7 > (l - p (W)) V | > Z - A)(W) by monotonicity of /, 
hence, /(7 + A)(W)) = +00 by definition of Iq. Otherwise, we have 7(7) < R, 
which means that the left term dominates the maximum. 

In both cases we get 

-?7/(7 + A)(W)) <hminfelog inf P(X u ' y eG). 

Now consider the unbounded case Iq = +00. Recall M defined by (1.27), 
and let (3q(G) := sup; >0 M(f(i)) ' ^ ne existence of which was claimed in the 
statement. If (3q(W) < (3q(G), we can choose l\ such that M J^^ > (3q(W) 
and set 7 := (3o(W)M(f(li)). Using the same arguments as in the bounded 
case, we deduce that 

-77/(7) < max/ lim inf e log inf F(X u ' y G G), -nf(h)\. 

Since / is increasing and l\ > 7, we obtain 

-77/(7) < lim inf e log inf Pfl^eG). 

e-»0 r<EW,y£K,fi€V,ueS r (s,n) 

In both the bounded and the unbounded case we have found 7 = (3q(W)5(G) 
such that the desired bound holds: we have 5(G) = M(R) + 1 in the bounded 
case and 5(G) = M(f(l\)) in the unbounded case. Furthermore, the choise 
W = {ro} corresponds to (3o(W) = (3 and yields 70 (G) = (35(G), in complete 
analogy to the situation of the upper bound. The uniformity is also proved 
in exactly the same way as already shown for the upper bound. □ 

2. Exit and entrance times of domains of attraction. We continue to 
study asymptotic properties of diffusions with weakly periodic drifts given 
by the SDE 

(2.1) dXl = b(^,Xt \ dt + VsdWt, t>0,X^ = x £M d . 
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In this section we shall work out the effects of weak periodicity of the drift on 
the asymptotic behavior of the exit times of its domains of attraction. This 
will be done under simple assumptions on the geometry associated to it. So 
we will have to specify some assumptions on the attraction and conservation 
properties of b. Essentially, we shall assume that W d is split into two domains 
of attraction, separated by a simple geometric boundary which is invariant 
in time. Apart from that, we shall assume that the drift is pointing inward 
sufficiently strongly so that the diffusions will not be able to leave compact 
sets in the small noise limit. Let us make these assumptions more precise. We 
recall that, according to the Kramers-Eyring law (see, e.g., [9]), the mean 
time a homogeneous diffusion of noise intensity e needs to leave a potential 
well of depth | is of the order exp | . Nature therefore imposes the time scales 
T £ with which we have to work. For simplicity, we measure these scales in 
energy units: with \i > 0, we associate the time scale T e = exp We assume 
as before that b satisfies the local Lipschitz conditions (1.14) and (1.15), and 
that the growth of the inward drift is sufficiently strong near infinity which 
is expressed by (1.16). 

The additional conditions concerning the geometry of b are specified in 
the following. 

Assumption 2.1. The two-dimensional ordinary differential equation 
(2.2) ip s (t) = b(s,ip s (t)), t>0, 

admits two stable equilibria x_ and x+ in R rf which do not depend on s > 0. 
Moreover, the domains of attraction defined by 

A±(s) = |y € R d :<Ps(t) = b(s,ip s (t)) 

(2.3) 

and v? s (0) = y imply lim tp s (t) = x± > 

t— »oo J 

are also independent of s > and denoted by A± . They are supposed to sat- 
isfy A- U A + = M. d , and cL4_ = dA + . We denote by x the common boundary 
(see Figure 1). 

The asymptotic results concerning the exit and entrance time should re- 
main true in a more general setting, where the stable equilibrium points and 
the domains of attraction depend on s > 0. We stick to Assumption 2.1 for 
reasons of technical and notational simplicity. 

The main subject of investigation in this section is given by the exit times 
of the domains of attraction A±, provided that the weakly time inhomoge- 
neous diffusion starts near the equilibrium points x±. By obvious symmetry 
reasons, we may restrict our attention to the case of an exit from A_. As 
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Fig. 1. Domains of attraction. 

we shall show, this exit time depends on the quasi-potential, that is, on the 
cost function taken on the set of all functions starting in the neighborhood 
of and exiting the domain of attraction through \- For this reason, we 
introduce the one-periodic energy function 

(2.4) e(s) : = inf V s (x_,y) <oo for s > 0, 

w yex ~ 

which is continuous on K + . In the gradient case b(t,x) = —X7 x U(t,x), this 
function coincides with twice the depth of the potential barrier to be over- 
come in order to exit from j4_, that is, the energy the diffusion needs to 
leave A_. Therefore, scales [i — corresponding to the Kramers-Eyring times 
T £ = exp(j) according to the chosen parametrization — at which we expect 
transitions between the domains of attraction must be comprised between 

jii* := inf e(t) and // := supe(t). 

These two constants are finite and are reached at least once per period since 
e(t) is continuous and periodic. Now fix a time scale parameter \x. This 
parameter serves as a threshold for the energy, and we expect to observe an 
exit from A_ at the first time t at which e(t) falls below (j,. For fj, g]/x*,/x*[, 
we therefore define (see Figure 2). 

(2.5) a u = inf{t>0:e(t) < n}, a M = inf{t > : e(t) < //}. 

The subtle difference between and may be important, but we shall 
rule it out for our considerations by making the following assumption. 

Assumption 2.2. The energy function e(t) is strictly monotonous be- 
tween its (discrete) extremes, and every local extremum is global. 
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Fig. 2. Definition o/a M and a^. 

Under this assumption, we have = a^. We are now in a position to 
state the main result of this section. Let g > be small enough such that 
the Euclidean ball B s (x+) C A + , and let us define the first entrance time 
into this ball by 

(2.6) T e = m£{t>0:Xf€B e (x+)}. 

This stopping time depends of course on e, but, for notational simplicity, we 
suppress this dependence. 

Theorem 2.3. Let f/, < e(0). There exist 7] > and ho > such that, for 
h < ho, 

limelog sup F y (T g (£{(a l ,-h)T £ ,(a u + h)T £ }) = n-e{a fl -h). 

yeB v (x.) 

Moreover, under Assumption 2.2, this convergence is uniform w.r.t. fx on 
compact subsets o/]/i*,e(0)[. 

Note that Assumption 2.2 implies the continuity of \i *— ► fx — e{a^ — h), and 
that ji — e{a a — h) < if h is small enough. The statement of the theorem 
may be paraphrased in the following way. It specifies time windows in which 
transitions between the domains of attraction will be observed with very 
high probability. In particular, if e(t) is strictly monotonous between its 
extremes, we prove that the entrance time into a neighborhood of x + will 
be located near a^T £ in the small noise limit. The assumption fi < e(0) is 
only a technical assumption in order to avoid instantaneous jumping of the 
diffusion to the other valley. It can always be achieved by simply starting 
the diffusion a little later. We could even assume that e(0) = //*, which then 
would yield uniform convergence on compact subsets of 

The rest of this section is devoted to the proof of this main result and is 
subdivided into separate subsections in which lower and upper bounds are 
established. 
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2.1. Lower bound for the exponential exit rate: diffusion exit. We have 
to establish upper and lower bounds on the transition time T e which both 
should be exceeded with an exponentially small probability that has to be 
determined exactly. It turns out that the probability of exceeding the upper 
bound (a« — h)T £ vanishes asymptotically to all exponential orders, so the 
exact large deviations rate is determined only by the probability W x (tq < 
(an — h)T £ ) of exceeding the lower bound. 

For a lower bound of the latter probability, as well as for an upper bound 
on P x (r e > (q m + h)T 6 ), one has to prove large deviations type upper bounds 
of the asymptotic distribution P x (r e > s(e)) for suitably chosen This 
can be expressed in terms of the problem of diffusion exit from a carefully 
chosen bounded domain. 

Recall that t q is the first entrance time of a small neighborhood B g (xj r ) 
of the equilibrium point x + . Consider for R, g> the bounded domain 

D = D(R,g) :=BM&)\B g (x + ), 

and let 

T D :=m{{t>0:X t (£D} 

be the first exit time of X from D. An exit from D means that either X 
enters B g (x^), that is, we have a transition to the other equilibrium, or 
X leaves -Br(O). But, as a consequence of our growth condition (1.16), the 
probability of the latter event does not contribute on the large deviations 
scale due to Proposition 1.4, as the following simple argument shows. 

Let s(e,fi) = sT e for some s > 0. Since td = T e Actr, where an is the time 
of the diffusion's first exit from Br(0), Proposition 1.4 provides constants 
Ri,£i > s.t. for R> R±, e < ei, 

F x .(r e > s(e,/i)) < F x {{t s > s(s,n)} n {a R > s(e, fi)}) + P x (a R < s(e,/x)) 
<Px(TD>s(e,fi))+C V 2 ^^ e~^ R l e for ||x|| < -. 

£ 2 

By the choice of s(e,fi) and T e = exp(^), the right-hand side term in the 
last sum is of the order g exp ^~^ R ; that is, it can be assumed to be ex- 
ponentially small of any exponential order required by choosing R suitably 
large. Obviously, this holds uniformly with respect to jjL on compact sets. 
This argument shows that the investigation of asymptotic properties of the 
laws of T e may be replaced by a study of similar properties of td, with an 
error that may be chosen arbitrarily small by increasing R. 

Similarly to the time homogeneous exit problem, we need a lemma which 
shows how to approximate the energy of a transition by the cost along 
particular trajectories which exit some neighborhood of D. This is of central 
importance to the estimation of the asymptotic law of tq. 
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Lemma 2.4. Let $ > and M a compact interval o/M+. Then there 
exist To > and 5 > w'i/i £/ie following property: 

For all x G D and s G M , we can find a continuous path C x ' s G Cor start- 
ing in Co' s = x and ending at some point of distance d(C^,D) > 5 away 
from D such that 

Iqt (C' s ) < e(s) + t? for all s G M. 

Proof. This proof extends arguments presented in Lemmas 5.7.18 and 
5.7.19 in [4]. 

Fix $ > 0, and let us decompose the domain D into three different ones. 
Fixing I > 0, define a domain Pi by Pi = {x G D:dist(x,x) < 0- We recall 
that x is the separation between A- and j4 + . Then we define two closed sets 
T>_ = (D \ Pi) n and D + = (D\Pi)f]A + . We shall construct appropriate 
paths from points y G D to points a positive distance away from D not 
exceeding the energy e(s) by more than ?? uniformly in s G M in four steps. 

Step 1 . Assume first that y G D- . For / > small enough, we construct 
S\ > 0, 5"' > and a path t/;^' defined on a time interval [0,r^ ,2/ ' ] with 
T s,y,i £ or a jj ^ g D_. s G M and along which we exit a ^-neighborhood 
of D_ at cost at most e(s) + |$. 

Step 1.1. In a first step we go from y to a small neighborhood L>i(x_) of 
x_, in time at most T{ < oo, without cost. 

We denote by ip*' y ' the trajectory starting at ip^' v ' (0) = y G D_ of 

^i(t) =&(s,<£i(t)), 

and reaching Bi(x~) at time a y,s ' 1 . Since D_ C A^ and due to Assump- 
tion 2.1, u^' 8 '' is finite. Moreover, since b is locally Lipschitz, stability of 
solutions with respect to initial conditions and smooth changes of vector 
fields implies that there exist open neighborhoods W y of y and W s of s and 
T?' y > 1 > such that, for all z G W y , u G W s , a?*' 1 < T*> y ' 1 . Recall that Z)_ 
is compact. Therefore, we may find a finite cover of D_ x M by such sets, 
and, consequently, T[ < oo such that, for all y G Z)_ and s G Af, <jj' y '' < T-f. 
Denote z 5 '^ = tp\' v ' (cr 8 ' 27 ' ). 

Step 1.2. In a second step, we go from a small neighborhood Bi(x_) of 
x_ to the equilibrium point x_, in time at most 1, at cost at most |. In fact, 
by the continuity of the cost function, for / small enough, s G M, there exists 
a continuous path ip^' °f time length a^' < 1 such that ^'^''(O) = z s,y ' 1 , 
ip s 2 ' v ' l (a s 2 ' v ' l )=x_ and / .wig* 1 ) < #/3. 
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Step 1.3. In a third step, we exit some ^-neighborhood of D_, starting 
from the equilibrium point X-, in time at most T3 < 00, at cost at most 
e(s) + § for s G M. 

By (2.4) and the continuity of the cost function for any s G M, there 
exists z s ^ D T| < 00, some neighborhood W s of s and for it G W s , 
we have 933 G Co<t" such that ^3(0) = x_, ^(a^) = z s , < T| and 

sup iz,M)<e(8) +m 

Use the compactness of M to find a finite cover of M by such neighborhoods, 
and, thus, some T3 < 00 such that all the statements hold with 03 < T3 for all 
s G M. Finally, remark that the exit point is at least a distance 5 = infj g j \zi\ 
away from the boundary of £)_, if Zi,i G J, are the exit points corresponding 
to the finite cover. 

In order to complete Step 1, we now define a path ijj^' 1 G C 0t s, v ,i by con- 
catenating tp s { y ' 1 1 tp2 V ' 1 and (/?§. This way, for small / > 0, we find S[ > such 
that, for all s G M, y G D_, we have r^'' < 5(, ip[^ l (T^ 1 ) = y,tpi' y ' l (r"' v ' 1 ) £ 
A- and 

< e(s) + ftf for all s G M, y G D_. 

At this point, we can encounter two cases. In the first case ijj s ^ v ' 1 exits a 
^-neighborhood of Br(0). In this case we continue with Step 4. In the 
second case, ij)^' 1 exits D_ into fli, and we continue with Step 2. 

Step 2. For I small enough, we start in y G j3i to construct 5^ > and 
a path V^'' defined on a time interval [0, t%' Vj1 ] with t^' 1 < S l 2 for all y G 
L>_ , s G M and along which we exit into the interior of D + at cost at most 

t 
3 • 

In fact, due to the continuity of the cost function (see Lemma 1.6), there 
exists / > small enough such that, for any s G M, y G /?/, there exists 
z S) yj in the interior of D + , such that, ^'^'(O) = y V^*^ 7 "!' 2 ^) = z s,y,i and 
r n W'') < tf/ 3 - We ma y take S 2 = !■ 

S'tep 3. We start in y G D + to construct <5 3 > 0, S3 > and a path V ; 3' 2/ ' Z 
defined on a time interval [0, r^' ] with t^' v ' < S 3 for all y G D_,s G M and 
along which we exit D+ into B g _ s i^ (&+) at no cost. 

Let 5 l 3 = q/2. Since D + is compact and contained in the domain of at- 
traction of x+, stability of the solutions of the differential equation ip{t) = 
b(s,(p(t)) with respect to the initial condition y G D + and the parameter s 
guarantees the existence of some time S3 > such that the entrance time 
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t^' v ' of B Q i 2 {x + ) by the solution starting in y is bounded by S3. There- 
fore, we may take ip^' 1 to be defined by this solution restricted to the time 
interval before its entrance into B g / 2 (x+). 

Step 4. For I > small enough, we start in x G TJ_ and construct Tq > 
0, 5 > and a path ^ s ' x defined on the time interval [0,To], exiting a 6- 
neighborhood of D at cost at most e(s) + 1? for all s 6 M. 

For / small enough, take Tq = S[ + S l 2 + S l 3 . We just have to concatenate 
paths constructed in the first three steps. Recall that ip*' x ' 1 passes through 
the equilibrium x_ due to Step 1. In case ip^ x ' exits a (^-neighborhood of 
Br{0), just let the path spend enough time in x_ without cost to obtain 
a path £ s ' x '* defined on [0, To], and take S = d[. In the other case, we con- 
catenate three paths constructed in Steps 1-3 to obtain a path defined on a 
subinterval of [0, To] depending on s,x,l and which exits a ^-neighborhood 
of D. Recall from Step 1 that this path also passes through x_. It remains 
to redefine the path by spending extra time at no cost in this equilibrium 
point to complete the proof. □ 

We now proceed to the estimation of uniform lower bounds for the asymp- 
totic law of td- The uniformity has to be understood in the sense of Re- 
mark 1.9. 

Proposition 2.5. Let K be a compact subset of D: 

(a) If e(s) > fi, then 

liminf elog inf ¥ x (td < sT e ) >fi- e(s), 

locally uniformly on {(s,fi) min(e(0), e(s)), < s < 1}. 

(b) If e(s) < fi, then 

lime log supP x (T£) > sT e ) = —00, 

locally uniformly on {(s,/x) :e(s) < fj, < e(0),0 < s < 1}. 

Proof. We choose a compact subset L of [0, 1] and a compact subset 
M of ]/i*, e(0)[, as well as some 1? > such that 

|e(s)-//|>i? V(s,/i)axM. 

Later on we shall assume that e(s) — [i is uniformly positive respectively 
negative in order to prove (a) respectively (b). 

In a first step, we apply Lemma 2.4 to approximate the energy function 
e(s) by the rate function along a particular path, uniformly w.r.t. s. For the 
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chosen it yields Tq> and 5 > 0, as well as continuous paths C x ' s indexed 
by x G D and s G [0, 1] ending a distance at least <5 away from D such that, 
for all x G D and s G [0, 1] , 

In a second step, we use the Markov property to estimate the probability 
of exiting D after time sT £ by a large product of exit probabilities after 
time intervals of length independent of e and \x. Since for e > 0, [i G M, the 
interval [0, sT £ ] becomes arbitrarily large as e — > 0, we introduce a splitting 
into intervals of length v>Tq independent of e and \x. For G No, let t k = 
tk(s, e, fx) := sT e — ku. Then we have, for k G No and x G D, 

F x (r D >t k )=E x (l {TD > tk} l TD > tk+1 ) 

= E x (1{ TD >t k+1 }E [1{ TD > tk } | Ft k+1 } ) 

< V x (t d > t k+ i) sup¥ y>t (t d > v). 

ydD 

Here Py jS denotes the law of X s,y , defined by the SDE 

dX;> v = b(j±±,X? v }dt + y/idWt, t>0, X s 'y = y€R d . 

On intervals [0, v\ it coincides with the law of the original process X on 
[s,s + v\. Denoting q^[s,e, fi) := sup y£D f y> t k (td > an iteration of the 
latter argument yields 

JV(e,A0 

(2.7) supP x .(r D >sT £ )< ' q k (a,e,n) 

X ^ K k=i 

whenever N(e,[i)i> < sT £ . For the further estimation of the q k , we apply 
some LDP to the product (2.7). This relies on the following idea. We choose 
N(e, jj) of the order eT e . Then the starting times t k appearing in the product 
belong to some neighborhood of sT e that, compared to T £ , shrinks to a 
point asymptotically. Consequently, the family of diffusions underlying the 
product is uniformly exponentially equivalent to the homogeneous diffusion 
governed by the drift b(s, •). 

This will be done in the following third step. For x G D, s G [0, 1], let 

be the open 5/2-neighborhood of the path chosen in the first step, and let 

\£(x):= |J *(z,s). 
se[o,i] 
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To apply our large deviations estimates in this situation, note first that 
conditions concerning td translate into constraints for the trajectories of 
X s as figuring in the preceding section: due to the definition of *ff(x,s), the 
choice v > Tq and Lemma 2.4, we know that, for y G D, k < N(e, fi), if X tk ' y 
belongs to ^f(x), then for sure X tk ' y exits D before time v. Keeping this in 
mind, we may apply Proposition 1.8 to the neighborhoods 

S'>°(e, n) = [sT £ - vN{e, fi),sT £ + v] 

of sT £ . Each of the intervals [tk,tk + v] is contained in S s ' (e,[i). As men- 
tioned before, N(e,fi) is chosen of the order eT £ , and this can be done uni- 
formly w.r.t. [i £ M. More precisely, we assume to have constants < c\ < c 2 
such that aeT £ < N(e,n) < c 2 eT £ . Then 

diamS' s ' (e, /i) 
hm sup — = 0, 

e ^°se[o,i] lAt GA/ Te 

and by the large deviations principle of Proposition 1.8, we obtain the lower 
bound 

>-sup/ s To (C^)>- e (s)-f. 

y£K ^ 

We stress that this bound is uniform w.r.t. s in the sense of Remark 1.9, 
so we can find £o > independent of s such that, for e < Eq, /i£ M and 
k<N(e,n), 

1 - q k (s, e, n) = inf F yt (r D < v) 

y£D 

> inf F y>tj (T D <u)>exp\--(e(s) + ^]). 



From this, we obtain 



N ( s ^) / f 1 / ,Q\ -V \ N(e,/j,) 



sup¥ x {r D >sT £ )< f{ q k ( 8l e, f i)<(l-exp{--(e(8) + t- 

= exp|iV(e, jt) log ^1 - expj ~ (e(s) + 7j) }) } = : m ( e > A*)- 

Since log(l — x) < — x for < x < 1, we have 

m(e,fj.) < exp|-cie expj ~ - ^(e{s) + 7^} }■ 

In the fourth and last step, we exploit this bound of m(e,fi) to obtain the 
claimed asymptotic bounds. 
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In order to prove (a), assume that \i < e(s) for (s,/x) G L x M. Then the 
inner exponential approaches on L x M. Using the inequality 1 — e~ x > 
xexp(— 1) on [0,1], we conclude that there exists e\ G (0,eo) such that, for 
all e < ei and (s, fi) G L X M, 



elog inf F x (t d < sT e ) > elog(l - m(e,n)) 

x&K 



> elog^eci exp(— 1) exp|- (fj, — e(s) 

= — e + e log ci + e log e + // — e(s) — 

> fi — e(s) — i9. 



_ t9 
~ 2 

2 



For (b), assume \x > e(s) on L x M. Then 
elog supP^T/) > sT e ) < elogm(e,/i) 



< -ci£exp|-^f/x-e(s) - 0|— ^-oo. n 

As a consequence of these large deviations type results on the asymptotic 
distribution of td and the remarks preceding the statement of Lemma 2.4 
and Proposition 2.5, we get the following asymptotics for the transition time 
of the diffusion. 

Proposition 2.6. Let x G . There exists h >0 such that 

(2.8) liminf elogP x (T B < (a„ - h)T £ ) > p- e(a M - h), 

(2.9) limelogP :r (r >(a M + /i)T £ ) = -oo, 

for h< ho. Moreover, these convergence statements hold uniformly w. r. t. x 
on compact subsets of D and w.r.t. fj, on compact subsets of e(0)[. 

Proof. As the estimation based on Proposition 1.4 at the beginning of 
the section shows, we may derive the required estimates for td instead of 
T e , if R is chosen large enough. 

Let M be a compact subset of ]/z*,e(0)[. Then < a M < 1 for fj, G M 
which yields the existence of ho > such that the compact set := {a^ — 
h:fi£ M} is contained in ]0, 1[ for h < ho- Moreover, we have e(s) > /i for 
< s < dp due to the assumptions on e, uniformly w.r.t. (s,p) G Lh x M by 
the continuity of e. Hence, by Proposition 2.5(a), 

liminf elog inf F x (td < sT £ ) > fi — e(s), 

e— >0 x&K 
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uniformly on x M for all h < Hq. By setting s = — h, we obtain the 
first asymptotic inequality. The second one follows in a completely analogu- 
ous way from Proposition 2.5(b) since a M = a M and e(a^ + h) < /i for small 
enough h. □ 

2.2. Upper bound for the exponential exit rate. Let us next derive upper 
bounds for the exponential exit rate which resemble the lower bounds just 
obtained. We need an extension of a result obtained by Freidlin and Wentzell 
(Lemma 5.4 in [18]). 

Lemma 2.7. Let K be a compact subset of A^ \{x_}. There exist Tq > 
and c > such that, for allT >Tq, s £ [0, 1] and for each <p 6 Cqt taking its 
values in K, we have 



Let t(s,x) be the first exit time of the path (j) sx from the domain K. Since 
A- is the domain of attraction of x_ and since K is a compact subset of 
A- \ we obtain t(s,x) < oo for all x G K. 

The function t(s,x) is upper semi-continuous with respect to the vari- 
ables s and x (due to the continuous dependence of <p StX on s and x). Hence, 
the maximal value T\ := sup s6 r 0jl i x£K r(a,x) is attained. 

Let To = T\ + 1, and consider all functions <p £ Cqt with values in K. This 
set of functions is closed with respect to the maximum norm. Since there 
is no solution of the ordinary differential equation in this set of functions, 
the functional /q Tq reaches a strictly positive minimum on this set which is 
uniform in s. Let us denote it by m. By the additivity of the functional Iq T , 
we obtain, for T >Tq and ip £ Cqt with values in K, 



the hitting time of a small neighborhood of the equilibrium point x + . First 
we shall consider upper bounds for the law of this time in some window 
of length (3T £ where (3 is sufficiently small. The important feature of the 
following statement is that (3 is independent of s, while the uniformity of 
the bound again has to be understood in the sense of Remark 1.9. 



I s QT {^)>c{T-T ). 



Proof. Let <j) S)X be the solution of the differential equation 
<j>s,x(t) = b(s, <j> s ,x(i)), (ps,x(0) = x£K. 




withc = ^. □ 

Ml 

Let us recall the subject of interest of this subsection: 
T e = M{t>0:XleB e (x+)}, 
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Proposition 2.8. For all ■& > 0, there exist (3 > 0, n > such that, for 
all s G [0,1], 

limsupelog sup F x (sT £ <t q <(s + (3)T e ) < // - e(s) + 0. 

JTiis bound holds locally uniformly w.r.t. \x G e(0)[ and uniformly w.r.t. 
a €[0,1]. 

Proof. Let M be a compact subset of ]/i*,e(0)[, and fix -# > 0. We first 
introduce some parameter dependent domains, the exit times of which will 
prove to be suitable for estimating the probability that r g is in a certain 
time window. 

For this purpose, we define for 5 > and s G [0, 1] an open domain 

D(6, s) := |y G R d : y) < ,u* + dist(y, A+) > 5 

and we let D = D(5) = U s g[o i] D(5, s). Then D is relatively compact in A_, 
dist(y,j4+) > <5 for all y G D(6), and a transition to a neighborhood of 
x + certainly requires an exit from D(5). The boundary of D(5) consists of 
two hyper surfaces, one of which carries an energy strictly greater than /i* 
and thus greater than e(s) for all s G [0, 1]. The minimal energy is therefore 
attained on the other component of distance 5 from A + which approaches 
X = dA_ as 5 — > 0. Thus, by the joint continuity of the quasi-potential, we 
can choose 5q > and n > such that, for 5 < 5o and s G [0, 1], 

e(s) = inf V s (x-,z) > inf V s (x-,z) 

■zex z£dD(6) 

(2.10) 

> inf inf V s (y,z)>e(s)--. 

y£B v (x-) zedD(6) 4 

Let td be the first exit time of X s from D. For s G [0, 1] and (3 > 0, we 
introduce a covering of the interval of interest [sT £ , (s + Z?)? 15 ] into N = 
N(/3,£,fi) intervals of fixed length v, that is, v is independent of e, /x, s and 
/?. We will have to assume that v is sufficiently large, which will be made 
precise later on. Thus, we have Nv > (3T £ , and we can and do assume that 
N < (5T e . For k G Z, jfe > -1, let 

be the starting points of these intervals. We consider t_i since we need 
some information about the past in order to ensure the diffusion to start in a 
neighborhood of the equilibrium x~. Then for x G B„(cc_), we get the desired 
estimation of probabilities of exit windows for t q by those with respect to 
Td" 

N 

W x (sT £ <r e <(s + I5)T S ) < < TD < tfc+i). 

fc=0 
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In a second step we will fix k > and estimate the probability of a first exit 
from D during each of the intervals [t k , t k +i] separately. Here the difficulty is 
that we do not have any information on the location at time t k . We therefore 
condition on whether or not X s has entered the neighborhood B„(x-) in 
the previous time interval. For that purpose, let 

o k : = inf{t > t k V : X £ t G B v (x-)}, k > -1. 

Then for k > 0, 

< r D < t k+1 ) < F x (t k <t d < i fe+ i|ajfe_i < t k ) 

(2.H) 

+ F x (t d AffH >t fc ). 

In the next step we shall estimate the second term on the right-hand side 
of (2.11). Let K = K(5,rj) = D(5)\B t) (x-). Then K is compact, and by the 
Markov property, we have 

Wx(t d A cr fc „! > t k ) < supP^t^j (r D A a x > u), 

y&K 

where P^j is as defined in the previous section. Now we wish to further 
estimate this exit probability using large deviations methods. The neighbor- 
hoods 

S s 'P{e, fi) = [sT £ - v, sT £ + vN((3, e, n)] 
of sT £ contain each interval [t k ,t k+ i], —l<k< N(/3,s, fx), and they satisfy 



limsup sup — < p. 

e^o ^eAf,se[o,i] T £ 



Hence, by the uniform LDP of Proposition 1.8, applied to the closed set 

$ K (8, r)) = {<P£ Co,v ■ <Pt G K(5, rj) for all t e [0, u]}, 
we obtain the upper bound 

limsupelog sup P lJ) t k _ 1 (td A a\ > v) 

£->0 y£K,fi€M,k<N 

(2.12) <limsupelog sup P y , t (X £ e <S> K (5,i])) 

<-inf inf Iq v (v), 

y€K ip£>5> K (S,r))-<o(P) ' 

where 7o(/3) = P5(^k(8,t])) is the "blowup-factor" induced by the diameter 
(5. Since 70 (/3) — ► as f3 — > 0, we can find (3o > such that, for f3 < Po, 

<MM 7o(/3) c<i>// 



2' 2 ' 
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which amounts to saying that, instead of blowing up the set of paths, we 

' <5 2 s 
-2' 2- 



consider the slightly enlarged domain K (|, 3). Thus, 



By Lemma 2.7, the latter expression, and therefore the r.h.s. of (2.12), ap- 
proaches — oo as v — ► co, uniformly w.r.t. s £ [0,1]. So the second term in 
the decomposition of P x (ifc < to < ifc+i) can be neglected since it becomes 
exponentially small of any desired order by choosing v suitably large. 

In the next and most difficult step, we treat the first term on the right- 
hand side of (2.11). It is given by the probability that, while X s is in B„(x-) 
at time cr/c-i, it exits within a time interval of length tk+i — <Jk~i < 2v. Hence, 
by the strong Markov property, 

^x(tk < TD < tfc+lkfe-1 < t k ) < SUp Py,t(TD <2u). 

tk-\<t<tk,y£B n {x-) 
Applying the uniform LDP to the closed set 

F D (5) := {ip £ Co^'.ipo G D(5),<p to $ D(5) for some to < 2v) 
yields the upper bound 

limsupelog sup Pj,^ (td < 2v) 

e^O y eB r) {x-),neM,teS s >l 3 {e,ijL) 

(2.13) 

<— inf inf I^uWi 

i/GB^(x_) ¥ , G F D (5)To(/3) ' 

where 70 (/3) = 2(36(Fd(S))- By the same reasoning as before, we can replace 
the blow-up of the path sets Fd(S) by an enlargement of the domain D(5). 
We find fa>0 such that, for < fa 

~ M x inf ,Jhv&)<- inf mi tt 2u (y) 

y£B v (x_) ^sFo^toOS) y£B v (x^)<p£F D {8/2) 

<— inf inf V s (y,z). 

y&B rl {x-)zedD(8/2) 

Now we apply (2.10) and recall the uniformity of the LDP w.r.t. s. We find 
£0 > such that we have, for e < eq, s £ [0,1], |i£M and (3 < (3i, 

elog sup Wy t t(T D <2u) 

yeB v (x-),teS«>P(e,n) 

(2.14) <- inf inf V s (y,z) + - 

y&B r ,(x_)z£dD(8/2) 4 

<-e(s) + -. 
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We finally summarize our findings. We conclude that there exists s\ > such 
that, for e < ei, fx G M and s G [0, 1], we have 

elog sup F x {sT £ <r e <(s + P)T £ ) 

x£B v (x-) 

<elog< V sup F x (t k < t d <t k+l \a k ^i <t k ) \ + - 

I k=0 xeB v (x-) J 4 



<elog<j/3T £ exp( — 



3 

+ - = elog/3 + ^-e(s) + -i? 



</i-e(s) 
This completes the proof. □ 

Remark 2.9. If we stay away from s = 0, in the statement of Propo- 
sition 2.8 the radius of the starting domain B rj (x-) can be chosen inde- 
pendently of the parameter It may then be brought into the following 
somewhat different form. 

Proposition 2.10. Let L and M be compact subsets o/]0, 1] respectively 
]/i*,e(0)[. Let n > be small enough such that B rj (x-) belongs to the domain 

{y G R d : V s {x_,y) < // for all s G L}. 

Then, for all $ > 0, there exists some (3 > such that we have 

limsupelog sup F{sT e < t q < (s + /3)T £ ) < \i - e(s) + 

£—>0 xeB v (x-) 

uniformly w.r.t. s G L and \i G M. 

Proof. To prove Proposition 2.10, one has to modify slightly the pre- 
ceding proof. Instead of just 77, one has to choose two different parameters: 
r/o for the definition of the starting domain D and some rj\ for the descrip- 
tion of the location of the diffusion at time t k , that is, for the definition of 
the stopping times a k . □ 

In the following proposition, we derive the upper bound for the asymp- 
totic law of transition times, corresponding to the lower bound obtained in 
Proposition 2.6. 

Proposition 2.11. Let fj, < e(0), and recall from (2.5) the definition 
= inf{i > 0:e(i) < fi}. There exist 7 > and ho > such that, for all 
h < ho, 

(2.15) limsupelog sup F x (t s < (a^ — h)T £ ) < fi — e(a fl — h). 

£^0 x&B- 1 {x-) 

This bound is uniform w.r.t. fi on compact subsets o/]/i», e(0)[. 
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Proof. Let M be a compact subset of ]/x*,e(0)[. To choose ho, we use 
our assumptions on the geometry of the energy function e. Recall Assump- 
tion 2.2 according to which e is strictly monotonous in the open intervals 
between the extrema ]//*,//*[. It implies that e is monotonically decreasing 
on the interval [a e (o),a^] for any /i G M. By choice of M, we further have 
a e / ) < inf^gM <V- Now choose /io such that inf^gjv/ cl^ — ho > a e r y Then we 
have, for h< ho, 

(2.16) inf a„-/t>0, 

(2.17) e(0)> sup e(a M -/i), 

(2.18) e(s) > e(a^ — /i) for all s <a^ — h. 

To see (2.18), note that, for < s < a e ( ) , by definition of a e ( ) , the inequality 
e(s) > e(0) > e(a^ — h) holds, while, for a e ( ) < s < a M — h by monotonicity, 
e(s) > e(a M - /i). 

Next fix h < ho. For /i G M, let Ao = Ao(/i) = 0, and Ai(/x) < inf AlG M(o/x — 
h)T £ . For iVGN*,we set A;0u) = Ai + (K - h)T £ — Ai), 2<i< N, thus 
splitting the time interval [0, (a M — /i)T e ] into the N intervals [A,(//), Aj + i (//)] , 
< i < N - 1. Then for 7 > 0, x G B 7 (x_), 

AT-l 

Px(r, < (a„ - /i)T e ) < ^ P x .(r, G [A^), A i+1 (/i)]), 

i=0 

which implies 

limsupelog sup P^r^ < (a M — /i)T £ ) 

e— »0 igB 7 (i_) 

< max limsupelog sup F x (t 8 G [Aj(/x), A i+1 (/x)]). 

0<«<7V-1 e ^o i6B 7 (i_) 

Fix $ > such that, for /i < ^OiA* £ Af, we have e(0) > e{a fjL — h) + $. This 
is guaranteed by (2.17). We shall show that 

limsupelog sup F x (r e G [Aj(/z), A m (/i)]) < /i - e(a M - /i) + 

uniformly in < i < A — 1 and \i G M. 

Let us treat the estimation of the first term separately from the others. In 
fact, by Proposition 2.8, setting s = 0, j3 = A\/T £ , we may choose A, eo > 
and 70 > such that, for Ai < AT £ , e < eo, 7 < 70, /U G M , the inequality 

elog sup P x (T e e[A o (ji),A 1 (p)])<ii-e(0)+6 

x£Bj(x — ) 
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holds. Now we use the inequality e(0) > e(a fl — h) +*&, valid for all [i G M. 
Hence, there exists A > 0, Eq > and 70 > such that, for Ai < AT £ , e < eo, 
7<7o,Ate Af, 

elog sup ¥ x (r g G [Ao(/x), Ai(jLt)]) <fi — e(a M — h). 

x£Bj(x-) 

Let us next estimate the contributions for the intervals [Aj(//), Aj+i(//)] with 
i > 1. We use Proposition 2.8, this time with s = Aj(/u)/T e , (3 = -^-j- swp„ eM a 
By the definition of an, we get e(s) > e(a«) for all s < a«. By (2.18), we have 
e(s) = e(Ai(/i)/T £ ) > e(a M - h). By Remark 2.9, 

lim sup elog sup P x (r e G [Aj(//), A i+ i (//)]) < // - e(a M - h) + 1?, 

e— >0 sG_B 7 (x_) 

uniformly w.r.t. 1 <i < N and /Li G M. Letting 1? tend to 0, which implies 
that N tends to infinity and Ai tends to zero, we obtain the desired upper 
bound for the exponential exit rate. □ 

3. Stochastic resonance. Given the results of the previous section on 
the asymptotics of exit times which are uniform in the scale parameter fj,, 
we are now in a position to reconsider the problem of finding a satisfactory 
probabilistic notion of stochastic resonance that does not suffer from the lack 
of robustness defect of physical notions, such as spectral power amplification. 
We continue to study the SDE 

dX £ t = b(^,Xfj dt + y/edWu t > 0, X £ = x G M. d , 

introduced before, thereby recalling that the drift term b satisfies the lo- 
cal Lipschitz conditions (1.15) and (1.14) in space and time, as well as the 
growth condition (1.16). Moreover, b is assumed to be one-periodic in time 
such that T £ is the period of the deterministic input of the randomly per- 
turbed dynamical system described by X £ . 

In typical applications, b = —V X U is given by the (spatial) gradient of 
some time periodic double- well potential U (see [13]). The potential pos- 
sesses at all times two local minima well separated by a barrier. The depth 
of the wells and the roles of being the deep and shallow one change peri- 
odically. The diffusion X £ then roughly describes the motion of a Brownian 
particle of intensity e in a double-well landscape. Its attempts to get close 
to the energetically most favorable deep position in the landscape makes it 
move along random trajectories which exhibit randomly periodic hopping 
between the wells. The average time the trajectories need to leave a poten- 
tial well of depth | being given by the Kramers-Eyring law T £ = exp(|) 
motivates our choice of time scales T £ = exp(^) and also our convention to 
measure time scales in energy units /i. 
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The problem of stochastic resonance consists of characterizing the optimal 
tuning of the noise, that is, the best relation between the noise amplitude e 
and the input period T e — or, in our units the energy parameter [i — of the 
deterministic system which makes the diffusion trajectories look as periodic 
as possible. Of course, the optimality criterion must be based upon a quality 
measure for periodicity in random trajectories. 

In this section we shall develop a measure of quality based on the transi- 
tion probabilities investigated in Section 2 and with respect to this measure 
for fixed small e (in the small noise limit e — > 0), exhibit a resonance energy 
/io(e) for which the diffusion trajectories follow the periodic forcing of the 
system at intensity e in an optimal way. We shall in fact study the prob- 
lem in a more general situation which includes the double-well potential 
gradient case as an important example, and draws its intuition from it. The 
deterministic system y s (t) = b(s, (p s (t)),t > 0, has to satisfy Assumption 2.1, 
that is, it possesses two well separated domains of attraction, the common 
boundary of which is time invariant. In the first subsection we shall describe 
the resonance interval, that is, the set of all parameter values \i for which 
in the small noise limit trivial behavior, that is, either constant or continu- 
ously jumping trajectories, are excluded. The second subsection shows that a 
quality measure of periodic tuning is given by the exponential rate at which 
the first transition from one domain of attraction to the other one happens 
within a fixed time window around a^T 6 . This quality measure is robust, 
as demonstrated in the last subsection: in the small noise limit the diffu- 
sion and its reduced model, a Markov chain jumping between the domains 
of attraction reduced to the equilibrium points, have the same resonance 
pattern. 

3.1. Resonance interval. According to Freidlin [5], quasi-periodic hop- 
ping behavior of the trajectories of our diffusion in the small noise limit of 
course requires that the energies required to leave the domains of attraction 
of the two equilibria switch their order periodically: if e± denotes the en- 
ergy needed to leave A±, then e + needs to be bigger than e_ during part 
of one period, and vice versa for the rest. We assume that e± both satisfy 
Assumption 2.2 and associate to each of these functions the transition time 

a ^( s ) = inf{t > s : e±(t) < fj,}. 

The time scales \x for which relevant behavior of the system is expected 
clearly belong to the intervals 

h = 

Our aim being the observation of periodicity, we have to make sure that 
the process can travel back and forth between the domains of attraction on 



inf ej(t),sup ei(t) 



*€{-+}• 




the time scales considered, but not instantaneously. So, on the one hand, in 
these scales it should not get stuck in one of the domains. On the other hand, 
they should not allow for chaotic behavior, that is, immediate re-bouncing 
after leaving a domain has to be avoided. 

To make these conditions mathematically precise, recall that transitions 
from Ai to A^i become possible as soon as the energy ei needed to exit 
from domain i falls below [i which represents the available energy. Not to 
get stuck in one of A±, we therefore have to guarantee 

a > maxinf eAt). 

i=± t>0 

To avoid immediate re-bouncing, we have to assure that the diffusion cannot 
leave A_i at the moment it reaches it, coming from A{. Suppose we consider 
the dynamics after time s > 0, and the diffusion is near i at that time. Its 
first transition to A_i occurs at time a^(s)T £ , where a„(s) is the first time 
in the original scale at which ej falls below fi after s. Provided e_j(a^(s)) is 
bigger than /t, it stays there for at least a little while. This is equivalent to 
saying that, for all s > 0, there exists 5 > such that on [a^(s), a^(s) + 5] we 
have fj, < e_j. Since by definition for t shortly after a^(s), we always have 
e i(t) < A 4 ; our condition may be paraphrased by the following: for all s > 0, 
there exists 5 > such that on [a^(s), a^(s) + 5] we have \i < maxj = ±ej. 
This, in turn, is more elegantly expressed by 

a < inf maxej(t). 
t>o i=± 

Our search for a set of scales fi for which the diffusion exhibits nontriv- 
ial transition behavior may be summarized in the following definition. The 
interval 

/* = 

is called resonance interval (see Figure 3). 

In this interval, for small e, we have to look for an optimal energy scale 
//(e) in the following subsection. See [9] and [8] for the definition of the 



max inf e,- (t) , inf max a (t) 

i=± t>o n ' t>o i=± v ; 
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corresponding interval in the one-dimensional case and in the case of two 
state Markov chains. In Freidlin's [5] terms, stochastic resonance in the sense 
of quasi-deterministic periodic motion is given if the parameter /i exceeds 
the lower boundary of our resonance interval. 

Let us consider the potential gradient case. Assume that b(t, x) = —X7 x U(t, x), 
t > 0, x G where U is some time periodic double-well potential with time 
invariant local minima x± and separatrix. Then A- and A+ represent the 
two wells of the potential, x the separatrix. The energy e± is, in fact, the 
energy some Brownian particle needs to cross %. Freidlin and Wentzell [6] 
give the link between this energy and the depth of the well. 

Lemma 3.1. If D±(t) = inf y&x U(t,y) — U(t,x±) denote the depths of 
the wells, then e±(t) = 2D±{t) for all t>0. 

This link is the origin of the name "quasipotential." The minimal energy 
e is reached by some path which intersects the level sets of the potential 
with orthogonal tangents. This path satisfies an equation of the form 

ip s = V x U(t,(p s ), s G (-oo,T), ip T G x- 

The resonance interval is given by Ir = ] maxj = ± inft>o 2Di(t), 
inf t > max i= ±i 2 A (t) [. 

3.2. Transition rates as quality measure. Let us now explain in detail 
our measure of quality designed to give a concept of optimal tuning which, 
as opposed to physical measures (see [13]), is robust for model reduction to 
Markov chains just retaining the jump dynamics between the equilibria of 
the diffusion. We shall use a notion that is based just on this rough transition 
mechanism. In fact, generalizing an approach for two state Markov chain 
models (see [8] ) , we measure the quality of tuning by computing for varying 
energy parameters fj, the probability that, starting in Xj, the diffusion is 
transferred to x_j within the time window [(a 1 — h)T £ , {a 1 + h)T e ] of width 
2hT £ . To find the stochastic resonance point for large T e (small e), we have 
to maximize this measure of quality in /x G Jr. The probability for transition 
within this window will be approximated by the estimates of the preceding 
section. Uniformity of convergence to the exponential rates will enable us to 
maximize in [i for fixed small e. 

Let us now make these ideas precise. To make sure that the transition 
window makes sense at least for small h, we have to suppose that > 0, i = 
± for /x G Ir. This is guaranteed if 



ej(0) > inf maxej(i), i = ±. 
t>0 i=± 
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If this is not granted from the beginning, it suffices to start the diffusion a 
little later. For g small enough so that B g (x±) C A±, we call 

(3.1) M(e,fi,g)=mm sup F x (rp € [(a* - h)T s , (a* + /i)T £ ]), 

l=± xeB e (a5i) 

transition probability for a time window of width h. Here 

ri = M{t>0:XteB e (xi)}. 

We are ready to state our main result on the asymptotic law of transition 
time windows. 

Theorem 3.2. Let M be a compact subset of Ir, ho > and g be given 
according to Theorem 2.3. Then for all h < ho, 

(3.2) limelog(l — M(e,fj,, g)) = maxju — ei(a l — h)}, 

s~^0 i=± P 

uniformly for fi € M. 

Proof. This is an obvious consequence of Theorem 2.3. □ 

It is clear that, for h small, the eventually existing global minimizer fiR,(h) 

of 

IrB max{/i - ei(a l - h)} 
i=± ^ 

is a good candidate for our resonance point. But it still depends on h. To 
get rid of this dependence, we shall consider the limit of nn(h) as h — > 0. 

Definition 3.3. Suppose that 

IrB /i i-> max{/i - a{a l - h)} 
i=± ^ 

possesses a global minimum HR(h). Suppose further that /j,r = lim^^o fJ-R(h) 
exists in Ir. We call fiR the stochastic resonance point of the diffusion {X s ) 
with time periodic drift b. 

Let us now illustrate this resonance notion in a situation in which the 
energy functions are related by a phase lag 4> G ]0, 1[, that is, e_(i) = e+(t + 
4>) for all t > 0. We shall show that in this case the stochastic resonance 
point exists if one of the energy functions, and thus both, has a unique 
point of maximal decrease on the interval where it is strictly decreasing (see 
Figure 4). 
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Theorem 3.4. Suppose that e_ is twice continuously differentiable and 
has its global maximum at t\, and its global minimum at ti, where t\ < £2- 
Suppose further that there is a unique point t\ < s < ti such that e_|u 1)S r 
is strictly concave, and e_ |] s ,i 2 [ * s strictly convex. Then fiR = e-(s) is the 
stochastic resonance point. 

Proof. As a consequence of the phase lag of the energy functions, 
max{/x - a{a l - h)} = {/J, - e_(a~ - h)}. 

Write a M = a~ and recall that on the interval of decrease of e„, a M = el (//). 
The derivative of [i i— > [i — e_ (a^ — h) has to vanish for the minimizer, which 
yields 

1 = e'_(a M -h)-a'^ = e'_(a M - h) • 1 . 

e_\a tl ) 

Our hypotheses concerning convexity and concavity of e_ essentially means 
that e"_ (s) = 0, and e'L \]ti,s[ < 0' e -l]s,t 2 [ > 0, which may be stated alterna- 
tively by saying that fj, \— > e'_ (a M ) has a local maximum at = s. Hence, for h 
small, there exists a unique point a^h) such that e'_(a M (/i) — h) = e'_(a At (/i)) 
and lim/i^o a /i(^) = s - ^° snow that a^h) corresponds to a minimum of the 
function (j, \— > [fj, — e_ (a M — /i)] , we take the second derivative of this function 
at a^(h), which is given by 

eLMh) - h)e'L(a^h)) - e'Lja^h) - h)eL(a^h)) 

But e / _(a /x (/i)),e'_(a At (/i)-/i) < 0, whereas e'!_{a^{h) - h) > 0,e'L(a^(h)) < 0. 
This clearly implies that a^(h) corresponds to a minimum of the func- 
tion. But by definition, as h — > 0, a^Qi) — > s. Therefore, finally, e_(s) is 
the stochastic resonance point. □ 

3.3. The robustness of stochastic resonance. In the small noise limit 
e — > 0, it seems reasonable to assume that the periodicity properties of the 
diffusion trajectories, caused by the periodic forcing the drift term exhibits, 
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are essentially captured by a simpler, reduced stochastic process: a continu- 
ous time Markov chain which just jumps between two states representing the 
equilibria in the two domains of attraction. Jump rates correspond to the 
transition mechanism of the diffusion. This is just the reduction idea ubiq- 
uitous in the physics literature, and explained, for example, in [11]. We shall 
now show that in the small noise limit both models, diffusion and Markov 
chain, produce the same resonance picture, if quality of periodic tuning is 
measured by transition rates. 

To describe the reduced model, let e± be the energy functions corre- 
sponding to transitions from to A± as before. Assume a phase locking 
of the two functions according to the previous subsection, that is, assume 
that e_(i) = e+(i + (p),t > 0, with phase shift <\> E ]0, 1[. So, let us consider a 
time-continuous Markov chain {Yf, t > 0} taking values in the state space 
S = { — ,+} with initial data Yq = —. Suppose the infinitesimal generator is 
given by 



G 



T £ J V T £ 
t \ ft 



where ip(t) = <p(t + (f>),t > 0, and ip is a 1-periodic function describing a rate 
which just produces the transition dynamics of the diffusion between the 
equilibria ±, that is, 

(3.3) <p(t) = exp{-^H i>0. 



Note that by choice of ip, 

(3.4) ip(t) = expj - | , t>0. 

Transition probabilities for the Markov chain thus defined are easily com- 
puted (see [8], Section 2). For example, the probability density of the first 
transition time a% is given by 



p(t) = ip(t) exp-l — / ip(s)ds'>, if i ■ 



(3.5) 



q(t) = ip(t + 4>) exp< — / ip(s + <j)) ds > , if i = +, 



o 



t > 0. Equation (3.5) can be used to obtain results on exponential rates of 
the transition times <Tj if starting from —i, i = ±. We summarize them and 
apply them to the following measure of quality of periodic tuning: 

(3.6) M(e,fi) = minPj(cr_j G [(of, - h)T £ , (a i + h)T £ \), e>0,^G/ R , 
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which is called transition probability for a time window of width h for the 
Markov chain. 

Here is the asymptotic result obtained from a slight modification of Theo- 
rems 3 and 4 of [8] which consists of allowing more general energy functions 
than the sinusoidal ones used there and requires just the same proof. 

Theorem 3.5. Let M be a compact subset of In, ho < sup(a~ 1 ,T/2 — 
a" 1 ). Then for < h < ho, 

(3.7) lim eln(l — jV(e,/i)) = max{/i — e_(o* — h)}, 

e— »0 i=± F 

uniformly for fi G M. 

It is clear from Theorem 3.5 that the reduced Markov chain Y e and the 
diffusion process X s have exactly the same resonance behavior. Of course, 
we may define the stochastic resonance point for Y £ , just as we did for X s . 
So the following final robustness result holds true. 

Theorem 3.6. The resonance points of (X £ ) with time periodic drift b 
and of (Y £ ) with exponential transition rate functions e± coincide. 

Acknowledgment. We are much indebted to an anonymous referee for 
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